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We calculate the structure of neutron star interiors comprising both the hadronic and the quark 
phases. For the hadronic sector we employ a microscopic equation of state involving nucleons and 
hyperons derived within the Brueckner-Hartree-Fock many-body theory with realistic two-body 
and three-body forces. For the description of quark matter, we use several different models, e.g. 
the MIT bag, the Nambu-Jona-Lasinio (NJL), the Color Dielectric (CDM), the Field Correlator 
method (FCM), and one based on the Dyson-Schwinger model (DSM). We find that a two solar 
mass hybrid star is possible only if the nucleonic EOS is stiff enough. 
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The possible appeai^ance of quark matter (QM) in the interior of massive neutron stars (NS) is 
one of the main issues in the physics of these compact objects. Calculations of NS structure, based 
on a microscopic nucleonic equation of state (EOS), indicate that for the heaviest NS, close to the 
maximum mass (about two solar masses), the central particle density reaches values larger than 
1/fm^. In this density range the nucleons start to loose their identity, and quark degrees of freedom 
play a role. 

The value of the maximum mass of a NS is probably one of the physical quantities that are 
most sensitive to the presence of QM, and the recent claim of discovery of a two solai^ mass NS 
[|I|] has stimulated the interest in this issue. Unfortunately, while the microscopic theory of the 
nucleonic EOS has reached a high degree of sophistication [^, the QM EOS is poorly known at 
zero temperature and at the high baryonic density appropriate for NS. One has, therefore, to rely 
on models of QM, which contain a high degree of arbitrariness. At present, the best one can do is 
to compare the predictions of different quark models and to estimate the uncertainty of the results 
for the NS matter as well as for the NS structure and mass. 

In this paper we will discuss a set of different quark models in combination with a definite 
baryonic EOS, which has been developed within the Brueckner-Hartree-Fock (BHF) many -body 
approach of nuclear- matter, comprising nucleons and also hyperons [^]. In particular, in section |] 
we review the derivation of the baryonic EOS in the BHF approach, whereas Section ^ is devoted 
to the relevant features of the hadron quai^k phase transition for the several QM EOS discussed. In 
section ^ we present the results regarding NS structure, combining the baryonic and QM EOS for 
beta-stable nuclear matter, and conclusions are drawn. 



1. EOS of nuclear matter within Brueckner theory 

The EOS constructed for the hadronic phase at T = is based on the non-relativistic Brueckner- 
Bethe-Goldstone (BBG) many -body theory [||], which is a linked cluster expansion of the energy 
per nucleon of nuclear matter, well convergent and accurate enough in the density range relevant 
for neutron stars. In this approach the essential ingredient is the two-body scattering matrix G, 
which, along with the single-particle potential U , satisfies the self-consistent equations 



U{k;p)= {kk'\G{p-e{k) + e{k'))\kk')a, (1.2) 



where V is the bare nucleon-nucleon (NN) interaction, p is the nucleon number density, CO is the 
starting energy, \kakh)Q{kaki,\ is the Pauli operator, e{k) =e{k;p) = ^k^ + U {k;p) is the single 
particle energy, and the subscript "a" indicates antisymmetiization of the matrix element. In the 
BHF approximation the energy per nucleon is 



^2 = :^ L {kk'\G{p-e{k) + e{k'))\kk')a. (1.4) 

^ k,k'<kF 



2 



Constraints on the quark matter equation of state from astrophysical observations 



G. F. Burgio 



The inclusion of nuclear thi^ee-body forces (TBF) is crucial in order to reproduce the correct satu- 
ration point of symmetric nuclear matter [Q ||]. The present theoretical status of microscopically 
derived TBF is quite mdimentary. Recent results have shown that both two-body and three- 
body forces should be based on the same theoretical footing and use the same microscopical pa- 
rameters in their construction. Results shown here were obtained with the Argonne Vig (V18) [||] 
or the Bonn B (BOB) potentials, and compared also with the widely used phenomenological 
Urbana-type (UIX) TBF [ |To| ] (in combination with the V18 potential). It should be stressed that in 
our approach the TBF is reduced to a density-dependent two-body force by averaging over the po- 
sition of the third particle, assuming that the probability of having two particles at a given distance 
is given by the two-body correlation function determined self-consistently. 



In the past years, the BHF approach has been extended with the inclusion of hyperons [11, 12, 



130 , which may appear at baryon density of about 2 to 3 times normal nuclear matter density. The 
inclusion of hyperons produces an EOS which turns out to be much softer than the purely nucleonic 
case, with dramatic consequences for the structure of the NS, i.e., the value of the maximum mass 
is smaller than the canonical value 1.44 Mq. The inclusion of further theoretical ingredients, such 
as hyperon-hyperon potentials [^] and/or three-body forces involving hyperons, could alter the 
baryonic EOS, but unfortunately they are essentially unknown. Another possibility that is able 
to produce larger maximum masses, is the appearance of a transition to QM inside the star. This 
scenario will be illustrated below. 

Starting from the EOS for symmetric and pure neutron matter, and assuming stellar matter 
composed of neutrons, protons, and leptons [Q], the EOS for the beta equilibrated matter can be 



obtained in the usual standard way [^, |14| |15|] : The Brueckner calculation yields the energy density 
of baryon/lepton matter as a function of the different partial densities, e (p„ ,Pp,Pe,p^). The various 
chemical potentials for the species / = n,p,e,ljL can then be computed straightforwardly. 



5£ 

dpi 



(1.5) 



and the equations for beta-equilibrium, 

Hi = biH,, - qiHe , (1-6) 
(bi and <7, denoting baryon number and charge of species /) and charge neutrality, 

J^Piqi = 0, (1.7) 

i 

allow one to determine the equilibrium composition {Pi{p)} at given baryon density p and finally 
the EOS, 

dp p dp 

In Fig. |l] we compare the EOS obtained in the BHF framework when only nucleons and leptons ai^e 
present (thick lines), and the corresponding ones with hyperons included (thin lines). Calculations 
have been performed with different choices of the NN potentials, i.e., the Bonn B, the Argonne V18, 
and the Nijmegen N93, all supplemented by a compatible microscopic TBF [S]. For completeness. 
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Figure 1: Pressure vs. the baryon number density of hadronic NS matter. Thick curves show results for 
purely nucleonic matter, whereas thin curves include hyperons. 



we also show results obtained with the Argonne V18 potential together with the phenomenological 
Urbana IX as TBF. 

We notice a strong dependence on both the chosen NN potential, and on the adopted TBF, the 
microscopic one being more repulsive than the phenomenological force. The presence of hyperons 
decreases strongly the pressure, and the resulting EOS turns out to be almost independent of the 
adopted NN potential, due to the interplay between the stiffness of the nucleonic EOS and the 
threshold density of hyperons [^. The softening of the EOS has serious consequences for the 
structure of NS, leading to a maximum mass of less than 1.4 solar masses [T3|], which is below 
the observed pulsar masses [16|. 



2. Quark Phase 

The properties of cold nuclear matter at large densities, i.e., its EOS and the location of the 
phase transition to deconfined QM, remain poorly known. The difficulty in performing first- 
principle calculations in such systems can be traced back to the complicated nonlinear and non- 
perturbative nature of quantum chromodynamics (QCD). Therefore one can presently only resort 
to more or less phenomenological models for describing QM. 

One of the most widely used approach is the MIT bag model [[l7|]. In the calculations shown 
below, we assumed massless u and d quarks, s quarks with a current mass of nis = 150 MeV, and 
either a fixed bag constant B = 90 MeVfm^^, or a density-dependent bag parameter. 



B(p) =B„ + (Bo-B~)exp 



(2.1) 



with B„ = 50 MeVfm , Bq = 400 MeVfm , and j8 = 0. 17. This approach has been proposed in 



[18], and it allows the symmetric nuclear- matter to be in the pure hadronic phase at low densities. 



and in the quark phase at large densities, while the transition density is taken as a parameter. We 
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find that the phase transition takes place at baryon density of the order of 2-3 Po, and that the mixed 
phase contains a small fraction of hyperons. 

The phase transition from the hadronic phase to the pure quark matter phase is usually per- 
formed through the Maxwell construction, or the more sophisticated Gibbs construction [[15]]. This 
is widely used for the study of neutron star structure, several numerical details can be found in 
many references [15|, and will not be repeated here. The important point to stress is that the fea- 
tures of the phase transition, such as its onset and the density interval over which it extends, depend 
crucially on the models used for describing the hadron and the quark phase. For example, we 
found that the phase transition constructed with the CDM [19] model is quite different from the 
one obtained using the MIT bag model. In the CDM, the onset of the coexistence region occurs 
at very low baryonic density. This implies a large difference in the structure of neutron stars. In 
fact, whereas stars built with the CDM have at most a mixed phase at low density and a pure quark 
core at higher density, the ones obtained with the MIT bag model contain a hadron phase, followed 
by a mixed phase and a pure quark interior. The scenario is again different within the Nambu- 
Jona-Lasinio model [pO|], where at most a mixed phase is present, but no pure quark phase. It has 
also been found that the phase transition from hadronic to QM occurs at high values of the baryon 
chemical potential when the Dyson-Schwinger model is used to describe the quaik phase. In some 
extreme cases, for particular choices of the parameters, no phase transition at all is possible. In 
fact, the Dyson-Schwinger EOS [ 21 ] is generally stiffer than the hadronic one, and the value of 
the transition density is high. We also found that with the DSM no phase transition exists if the 
hadronic phase contains hyperons, just like the phase transition with NJL model. 

It is worthwhile to mention that neutron stars observations can help to determine free parame- 
ters, which could be present in some quark matter models. This is the case of the QM EOS based 
on the Field Correlator method [22], which depends crucially on the value of the gluon condensate 
G2- It turns out that using value of G2 — 0.006 — 0.007 GeV*, which gives a critical temperature 
Tc ~ 170 MeV, produces maximum masses which are only marginally consistent with the obser- 
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vational limit, while larger masses are possible if the value of the gluon condensate is increased. 
Also in this case, the phase transition only takes place if no hyperons are present in the hadronic 
phase. 



3. Neutron star structure 



We assume that a NS is a spherically symmetric distribution of mass in hydrostatic equilibrium. 
The equilibrium configurations are obtained by solving the Tolman-Oppenheimer-VoUcoff (TOV) 
equations [0] for the pressure P and the enclosed mass m, 

dP _ Gme (1 +P/e) {\+4nr^P/}n) 
dr \—2Gm/r ' 

dm 9 

— = AK?e , (3.2) 
dr 

being G the gravitational constant. Starting with a central mass density e{r = 0) = Ec, we integrate 
out until the density on the surface equals the one of iron. This gives the stellar radius R and the 
gravitational mass is then 

rR 

MG = m{R)=47i dr?e{r). (3.3) 

JO 

We have used as input the EOS discussed above, and the results are plotted in Figs. § and |3[ where 
we display the gravitational mass Mq (in units of the solar mass Mq = 2 X lO^^g) as a function of 
the radius R and central baryon density . 

Calculations displayed in Fig. ^ for neutron star matter are obtained in the BHF theoretical 
framework with the VI 8 nucleon-nucleon potential in combination with i) the CDM (solid black 
curve), and ii) the MIT bag model (blue curves) for quark matter. Due to the use of the Maxwell 



construction, the curves are not continuous [|15||: for very small central densities (large radii, small 
masses) the stars are purely hadronic. We observe that the values of the maximum mass depend 
only slightly on the EOS chosen for describing quark matter, and lie between 1.5 and 1.6 solar 
masses. A clear difference between the two models exists as far as the radius is concerned. Hybrid 
stars built with the CDM are characterized by a larger radius and a smaller central density, whereas 
hybrid stars constructed with the MIT bag model are more compact. Further calculations have been 
performed using the Paris potential as nucleon-nucleon interaction for the hadronic phase, and the 
NIL model for the quark phase (red curve). In this case the onset of pure quark matter leads to an 
instability, as well as in the FCM (green curve). In both cases the phase transition takes place only 
if no hyperons are present in the hadronic phase. Unfortunately, for all the cases discussed above, 
the value of the maximum mass lies below the mass observed for the pulsar PSR J 16 14-2230, 
Mq = 1.97 ±0.04 Such a high value puts severe constraints on the EOS, and in particular it 
demands an additional repulsion in the QM EOS. 

In Fig. ^ we show results obtained in the BHF framework using the Bonn B nucleon-nucleon 
potential, which produces the stiffest EOS, as shown in Fig. |]. This yields a very high maximum 
NS mass, « 2.50 Mq, with only nucleons (solid black curve), and 1.37 Mr.-, including hyperons. 
Using the DSM for the quai^k phase, we found that the maximum mass of hybrid stars is only 
a little lower than 2.5 with a = 0.5, and decreases to about 2 Mq with a = 2, being a a 



6 



Constraints on the quark matter equation of state from astrophysical observations 



G. F. Burgio 




Figure 3: Gravitational NS mass vs. the radius (right panel) and the central baryon density (left panel) for 
different EOS employing the BOB hadronic model. 



model parameter which controls the rate of approaching asymptotic freedom. With increasing a 
we can obtain a smooth change from the pure hadronic NS to the results with the MIT bag model. 
Moreover, no phase transition can occur and no hybrid star can exist if hyperons are introduced. If 
hyperons are excluded, the phase transition from nucleon matter to QM takes place at rather large 
baryon density, and the onset of the phase transition is determined in this case by the parameter a. 

The possible effects of the hadron-quark phase transition are very different with the MIT bag 
model and the DSM: In the case of the MIT model, the phase transition begins at very low baryon 
density and thus effectively impedes the appearance of hyperons. Consequently the resulting max- 
imum mass of the MIT hybrid star is 1 .5 Mq, lower than the value of the nucleonic star, but higher 
than that of the hyperon star given before. A clear difference between the two models exists as far 
as the radius is concerned. Hybrid stars built with the DSM are characterized by a larger radius 
and a smaller central density, whereas hybrid stars constructed with the MIT bag model are more 
compact. 

In conclusion, a hybrid star with 2 Mq is only allowed if the nucleonic EOS is stiff enough, 
and the hadron-quark phase transition takes place without hyperons in the hadronic phase. 
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